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Abstract 

We calculate the one-loop contribution to the bino annihilation rate due to the process 
BB Z* , which vanishes at tree level. 



It is very likely that some extension of the standard electroweak model is needed. In fact, 
such an extension is probably dictated by cosmology in order to produce a baryon asymmetry 
and to provide a suitable dark matter candidate. The former typically relies on grand 
unification (though simpler extensions also work |jl|), while the best motivated extension 
which provides a cold dark matter candidate is the minimal supersymmetric standard model 
(MSSM). Because of an unbroken discrete symmetry, R-parity, the lightest supersymmetric 
partner (LSP) is expected to be stable. This makes the LSP a natural candidate for the 
role of the non-baryonic dark matter. In previous studies, the lightest neutralino, a linear 
combination of the supersymmetric partners of the neutral SU{2) gauge boson (neutral 
wino), U{1) hypercharge gauge boson (bino), and the two neutral components of the Higgs 
doublets (higgsinos) has been suggested as the most likely candidate for the LSP In a 
large portion of parameter space, the least massive eigenstate is a nearly pure bino 0. 

Neutralino annihilation rates have been calculated previously at tree level, under a variety 
of assumptions and to various degrees of accuracy p|, ||, ||, H, 0, ||, ||, [111. For neutralinos 



which are nearly pure bino, the rate is dominated by creation of fermion pairs via the 
exchange of a sfermion, though many of the calculations have included annihilation channels 
into Higgs scalars as well. Since binos are SU{2) x f/(l) singlets, they do not couple directly 
to the Z. Thus, s-channel annihilations through Z exchange which are very important for 



massive neutrino annihilations or higgsino co-annihilations |T^, |Tl|, do not occur for pure 
binos. Of course in reality the LSP is never a "pure" bino, but rather always contains 
some small admixture of higgsinos as well. Annihilations through Z exchange due to this 
admixture make a small contribution to the cross-section. In this paper, we consider a 
correction to the bino annihilation rate from the process BB Z*, via the loop diagrams 
shown in fig. (1). Generically, the corrections are found to be much larger than the corrections 
due to the higgsino admixture. However, their contribution to the total annihilation cross- 
section remains small, on the order of a few percent or less over the physically allowed range of 
parameters. Recently one- loop corrections to the neutralino masses have been computed [|T2| , 



and found to be very small (less than one per cent) in the parameter range of interest here. 
These corrections to the masses would make similarly small corrections to the annihilation 
cross section, and so are significantly smaller than those we compute here. 

In general, the lightest neutralino mass eigenstate is a linear combination of the gauginos 
and B, and the higgsinos and ifj- 

x'^ = aW^ + (3B + -fH^ + 6H^ . (1) 
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In the pure bino region {(3 ~ 1), m~o ~ Mi, where Mi is the supersymmetry breaking U{1) 
gaugino mass. For our purposes here, the supersymmetric parameter space can be described 
by one additional parameter, e, the supersymmetric Higgs mixing mass (often denoted by 
— /i). In the bino region, results are very insensitive to the ratio of Higgs vevs, tan/5 (tan/5 
should not be confused with the j3 in eq. (1)). In what follows we have chosen tan/5 = 2 
and assumed a top quark mass of 160 GeV. We simplify by using GUT boundary conditions 
to relate M2, the SU{2) gaugino mass, to Mi via Mi = |M2tan^6'„. From here on, we 
will assume that e ^ M2 so that we are working in the nearly pure bino region and write 
m~o = m~. 

We begin by calculating the annihilation rate of almost pure binos to fermion anti-fermion 
pairs FF via (1) t-channel sfermion exchange, (2) BB ^ Z* ^ FF via the loop diagrams 
of fig. (1), and (3) BB ^ Z* ^ FF via the small but non-zero admixture of and H2 in 
X° (see eq. (1)). 

The relevant interaction terms are 
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cos6'„ 
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where P± = (1 ± 75)/2, 1/, = 2Qf, Yj, = 2iQf - T- 



-Ti/2 + Qf sin' 9^, and 
3 — +1/2 for up-type quarks 
and neutrinos, and —1/2 for down- type quarks and charged leptons. The sum runs over all 
fermions. We neglect the effect of sfermion mixing, which is treated in detail in ref. H]. 



= —Tl/2. The fermion charge is denoted by Q/, and 
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To derive a thermally averaged cross section, we make use of the technique of ref. [|13 
We expand {crvrei) in a Taylor expansion in powers of x = T/m 



(3) 



The coefficients a and b are given by 



&F+ -3 + 



Snip 



(4) 



where the subscript F specifies the final state fermions, ap and bp are computed from the 
expansion of the matrix element squared in powers of p, the incoming bino momentum, and 
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vp — {i — rrip / m'~y/'^ is a factor from the phase space integrals. In terms of the squared 
reduced transition matrix element dp can be written simply as, 
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327rm2~ 
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\Tp{p = 0)\' 



(5) 



Denoting the reduced transition amplitudes for the processes (1), (2) and (3) listed above 
by %fex, Tioop and respectively, we have 
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and where 



/ 



(9) 



(10) 



/ B^ ' B ' f ' 
Here C^, Yp^, Yp^ are the quantum numbers of the final state fermion, and nip and are 
the masses of the final state fermion and its supersymmetric partner. These expressions for 
the amplitudes do not include spin dependence, but we note that they will vanish if either 
the binos or fermions are spin parallel. In the other cases the amplitudes are antisymmetric 
under exchange of spins in the incoming or outgoing states. So their relative signs remain 
the same allowing us to directly sum and average over spin states giving eq. (5). To compute 
bp, where we need to consider p 7^ 0, this simple spin dependence no longer holds and we 
use standard trace techniques to calculate the squared transition matrix element. 

To compare the relative effect of the loop and the higgsino mixing diagrams on the 
dominant sfermion exchange process, we need to find the ap and bp coefficients for 1) the 
sfermion exchange alone, 2) the sfermion exchange plus the loop diagram and 3) the sfermion 
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exchange plus the higgsino mixing diagram. The dp for each combination can be read off 
from eq. (5) by redefining Tp to contain different combinations of diagrams, but the bp 
will need to be separately computed in each case. The expressions are too complicated to 
reproduce here. 

We now discuss the effect of this loop correction on the rehc density and detectabihty 
of binos. In fig. (2), we show a contour plot of ft^h'^ as a function of mg and the common 
sfermion mass mj — m^, assuming that binos annihilate to fermions only through the 
sfermion exchange process (1). is the mass density of the binos in units of the critical 
density, and h is the Hubble parameter in units of 100 km s~^ Mpc~^. The region where nip > 
mg is excluded, since here the LSP would be a sfermion rather than the bino. Notice the 
effect of the top threshold at — 160 GeV. Because the cross-section is p-wave suppressed, 
we have ap oc {mp/m-), and so as we cross below = mtop the relic density rises 
sharply. Notice also that Qg/i^ < 1 requires <, 500 GeV, while Jlg/i^ < 1/4 requires 
nip ^ 350 GeV. Keep in mind that these are approximate values for Q^h^ as only the 
dominant part of the cross-section was used here, though the corrections are very small. 

In fig. (3), we compare the effect of including the loop diagrams (2) on the relic density of 
binos fl^h^ to the effect of including the higgsino mixing diagram (3). As e gets larger, the 
admixture of higgsino in x° gets smaller, and so does the prefactor (7^ — S^) in eq. (8) (and 
in the corresponding p ^ expression). In the exactly pure bino limit, corresponding to 
£ — > 00 at fixed Ml, the process x^X^ — ^ ^ is forbidden, as mentioned above. For each point 
in the — rrip plane we find the value of e at which the magnitude of the loop correction 
is equal to that of the higgsino mixing correction. We then plot several contours of these e 
values. For each s, it is the region to the left of the corresponding contour where the loop 
effect is greater. In the parameter range of interest, we find that the loop correction to fl^h"^ 
vanishes near the line — mg, which is why e tends toward large values there. For the 
moderate values of m~ represented in fig. (3), we see that the effect of the loop is larger 
than the effect of higgsino mixing for large ranges in in the pure bino regime. Again, 
we exclude the region where > m^. We also exclude the region below the top threshold 
where both diagrams have a negligible effect on Vt^h?. 

Although we have found that for a large portion of the parameter space the loop correction 
is larger than the commonly used correction due to a higgsino admixture, the net effect on 
the annihilation cross-section and ultimately on VL^^h? is small. For all values of and M2 
consistent with O^/i^ < 1, the effect is at most a few percent, and in the allowed region, the 
effect predominantly decreases the relic density. Below the top threshold, the effect of the 
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loop diagrams drops by a factor of ~ 300 due to the accumulation of a number of unrelated 
factors. (This is also true for the higgsino mixing diagram.) Of course in the very thin slice 
of parameter space near = the annihilation cross-section will be greatly enhanced, 

and the effect of the loop diagrams will be significant. Because this effect is limited to a 
small region of parameter space we do not consider it further. We note only that in such 



regions the expansion of the thermal average should be done as explained in ref. [T^ • 

In addition to the effect of the loop correction on the relic density, it is also of interest to 
compute the effect on the zero temperature annihilation cross-section (p = 0) which governs 
annihilations in the galactic halo. Indeed, galactic halo LSP annihilation products have been 
frequently discussed as a potential signatures for dark matter [Q. 

In fig. (4), we show the effect of the loop diagrams on the zero-temperature cross-section, 
((Tfrei), at e = oo. The contours show the percentage change in a due to the inclusion 
of the loop correction diagrams (2). The effect is larger than on the relic density, but in 
the physically allowed region the effect is still <, 11%. We again see the top threshold at 
mg = 160 GeV. Notice that below the top threshold, the sign of the effect changes. At zero 
incoming bino momentum, the coupling of the Z to the final state fermions is proportional to 
Tg^, as seen in eq. (9). Above the top threshold, zero temperature annihilation is primarily to 
tops, while below the top threshold it is primarily to 6's and r's. Thus '7ioop(p = 0) switches 
sign, while %ie^{p = 0) does not. 

Finally in fig. (5), we compare the effect of including the loop diagrams (2) on the zero- 
temperature cross-section a to the effect of including the higgsino mixing diagram (3) . Calling 
the change in a due to the loop diagrams AioopCt and the change in a due to the higgsino 
mixture AhmO, we plot Aioopfl/AhmO as a function of e and using = 500 GeV. As e 
becomes large, the effect of the loop correction greatly exceeds the effect of higgsino mixing 
except very near to the value of mg at which the loop's effect on the cross section vanishes. 

To summarize, we have calculated the effect of the loop diagram of fig. (1) on the annihi- 
lation cross-section for binos. We find that away from the Z-pole, the effect of the additional 
diagram is to change the relic density of binos Q-^h"^ by at most a few percent. This is true 
over the physically allowed range of common sfermion mass and bino mass mg, though 
it is larger than contributions from the direct annihilation to Z's via the small but non-zero 
admixture of higgsinos in typically included in this type of calculation for a wide range 
in values of the higgsino mixing mass e. Near the Z-pole, we expect the effect of loop cor- 
rections to be considerably larger. The effect of the loop diagram on the zero-temperature 
cross-section is larger and may be as large as 11%. 
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Figure Captions 



Fig. 1) Feynman diagrams for loop corrections to BB — > Z*. 

Fig. 2) Relic density fl^h^ of binos, assuming bino annihilation only due to sfermion 
exchange. We take tan (3 — 2, mtop = 160 GeV, and a common sfermion mass 
rrip — mj. 

Fig. 3) A comparison of the effects of the loop diagrams and the higgsino mixing 
diagram on the relic density of binos ^^h?-, for several values of e. In the 
area to the left of each curve labelled by e, the loop correction has a greater 
effect on Qg/i^ than the higgsino mixing correction. 

Fig. 4) The effect of the loop diagrams on the zero-temperature cross-section. Con- 
tours label the percentage change in a due to the inclusion of the loop cor- 
rection. 

Fig 5) A comparison of the effects of the loop diagrams and the higgsino mixing 
diagram on the zero-temperature cross-section. 
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Fig 3 
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